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^ ■ Abstract 

\o 

I A method is given for calculation of a distribution of small particles, embedded 

' in a medium, so that the resulting medium would have a desired radiation pattern 

■ for the plane wave scattering by this medium. 

o 



"o.: 1 Introduction 



Let Dq C M be a bounded domain, ko be the wave number in Dq and k < ko he the 
wave number in Dq = M.^ \ Dq. We assume that Dq is a homogeneous medium, so that 
kQ is a constant. This assumption can be weakened: we may assume that kQ = kQ{x) is a 
known function. Let Dm be a particle, dm be its diameter, 1 < m < M, M is the number 
^ ' of small particles, a = maxi<m<M % is an estimate for the radius of a small particle. 

d • We assume that /cqO <^ 1, and then /ca <^ 1, i.e., particles are small, that d ^ a, where 

d = mmj-^mdist{Dm, Dj) and that the particles are acoustically soft, i.e., u = on Sm, 
the boundary of Dm- Denote Um=iDm '■= U , \ U = V , dU is the boundary of U, Cm 
is the electrical capacitance of a perfect conductor with the shape of Dm- 
The Inverse Problem is: 

Can one distribute many small particles in Dq so that a plane wave uq := e*'^"'^, 
(where a is a given unit vector, a ^ S'^ , S'^ is a unit sphere), scattered by Dq, would 
produce a desired radiation pattern (scattering amplitude) A{a',a)? 
The acoustic pressure u solves the problem (l)-(3): 



(V^ + k^ - q{x))u = in y, q{x) 




kl in Do 



in D'q, 
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M = on dU, 



(2) 



u = e''""''' + A(a' a) hoi-], r := Ixl oo, a' = -, (3) 

r \r J r 

and the coefficient A{a' , a) is called the scattering amplitude or radiation pattern. 
Let G = G{x, y) solve the problem 

[V^ + - ^(x)]^ = - 2/) in M^ r - ikG^ = o{l) r ^ oo. (4) 
This G exists, is unique, and solves the integral equation 

I" ^ik\x-y\ 

G{x,y) = g{x,y) - {k"^ - kl) g{x, z)G{z,y) dz, g = g{x,y) := — -. (5) 

Jdo 47r|x — y\ 

For brevity, we drop the fc-dependence in g, G and other functions below. One may 
consider G{x,y) known since Dq and k^ are known. 
Let us look for the solution to (jl))-© of the form 

M „ 

u = Uq{x , a) + ''^^ / G{x,t)am{t) dt, (6) 

m=l ^ 

where am{s) are unknown functions and Uq{x, a) is the scattering solution corresponding 
to the potential q{x) in the absence of small bodies, i.e., in the whole space. For arbitrary 
(Tm the right-hand side of (0) solves equation (|T)) (because G solves (jU), satisfies the 
radiation condition Q, and 

1 [ 

A{a',a) = Ag{a',a) + — Uo{s, -a')am{s) ds, (7) 

m=l '^•5'™ 

where the scattering solution Uo{s,a) can be defined by the formula ([H p. 232]): 

G(x, s) = Uo{s, a) + I - I , r := |x| — s> oo, a = — . (8) 

Anr \r J r 

Formula (jH)) was proved in ||3j p. 46], where it was shown that Uq{x, a) solves the problem 



N'^ + k'^ -q{x)]Uo{x,a) = OmR\ f/o = e^'^"'^ + a)— + o ( - ) , (9) 

r \r J 

where r := \x\ — >■ oo along the direction a'. The function Uq can also be considered 
known, and G can be considered as Green's function for the Schrodinger equation 0. 
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The right-hand side of (jH)) solves problem (H))-© if and only if am are such that boundary 
condition (j2)) is satisfied: 

-[/o(s,a)-^/ G{s,t)am{t)dt= G{s,t)aj{t) dt, s e Sj, I < j < M. 

(10) 

So far the smallness of the particles was not used. If fca <^ 1, then equation (jlOp can be 
simplified: 

-[/o(sj, ot)-^ G{sj, tm)Qm = / 9o{sj, t)aj dt, I < j < M, Qm ■= dm dt, 

(11) 

where Sj G Sj is any point in Sj, and Qq = ^^|^_^| . In equation (11) we have used two 
approximations. The first one is 

G{s,t)^g{s,t), s,teS,. (12) 

This approximation is justified by © when k\x — y\ ^1, because the term g{x, y) is the 
main term on the right-hand side of © as x — > y, and \sj — t\ <^ 1. The second one is 
g ^ Qq if k\s — t\ <^ 1, and the error of this approximation is 0{ka), so this approximation 
is also justified because ka <^ 1. 

Equation (|TT| is equation for the charge distribution aj on the surface Sj of a perfect 
conductor charged to the potential —Uo{sj, a) — J2m^j ^i^j^ tm)Qm- Therefore the total 
charge Qj := ajdt on the surface Sj can be calculated by the formula 



Qj = I ~Uo{sj,a) - ^G{sj,tm)Qni I 



1<3<M, (13) 



where Cj is the electrical capacitance of the perfect conductor Dj. In |T, p. 385] formulas 
for calculation of Cj with arbitrary desired accuracy are given. Equation ()13|) is a linear 
algebraic system for finding unknown Qj, 1 < j < M. 

Consider the limiting case M — > oo of the distribution of particles in Dq. Define C{y) 
as follows: 

/ C{y) dy = \im V C^, (14) 

D^CD 

where D G Dq is an arbitrary subdomain of Dq. The above definition means that C{y) 
is the limiting density of the capacitances of the small particles in D. Formula (jU)) shows 
that the self-consistent field Ue can be defined as 

, , \Uo{x,a) + Y,m=iG{x,t mini<rn<M\x -tm\'> a, 

Ue{x,a) = <. ~ ~ (15) 

[Uo{x,a) + Y.m^jG[x,tm)Qm, \x-tj\r~u a. 
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Therefore, defining we neglect the influence of any fixed single small particle on the 
field. This is justified when M ^ oo. Using (fT^ and (fT^ one gets 

Ue{x, a) = Uo{x, a) - ^ G{x, tm)CmUe{tm, «), 

m 

and in the limit M —>■ oo one obtains the equation: 

Ue{x,a) = Uo{x,a) - G{x,y)C{y)ue{y,a) dy, (16) 

where C{y) is defined in fjl4j) . Equation (fTBj) is equivalent to the Schrodinger scattering 
problem: 

[V^ + P - q{x) - C{x)]ue = in R^, (17) 
where q{x) is known, the function Ue has the following asymptotics 

= e^^"-^ + /l(a',a)— + o(- 1 , r := |x| ^ cx), a' = -, (18) 
r \r J r 

and v4(a', a) is the scattering amplitude at a fixed A; > 0. 

Therefore the Inverse Problem, stated above, is reduced to inverse scattering problem 
of finding the potential q{x) + C{x) from the knowledge of the corresponding fixed-energy 
scattering amplitude A{a', a). 

This problem was solved by the author (see Pfl Chapter 5] and references therein) . In 
Section 121 we outline the author's algorithm for solving this inverse scattering problem. 

If the potential q{x) + C{x) is found and q{x) is known, then C{x) is found and one 
knows the density of the particle distribution in D which produces the desired radiation 
pattern A{a',a). Assuming that the particles are identical, one has Cm = C, where C is 
the electrical capacitance of one particle, and C{x) = N{x)C, where N{x) is the density 
of particles, that is, the number of particles per unit volume around point x. See j2] for 
the theory of wave scattering by small bodies. 



2 Solution to inverse scattering problem 

We follow p. 264] and take k = 1 without loss of generality. Given A{a', a) one finds 
A£{a) := J^2 A[a', a)Y£{a') da', where := 1^™, are the normalized spherical harmonics, 
so that 

oo oo oo £ 

e=o e=o 1=0 m=-£ 

If the data A(a',a) are exact, then (jU p. 262]) 

^\Ma)\<o(^{^f\ \Y,m<^£y^ Vr>0. , e M. (20) 
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where feo > is the radius of the smallest ball containing the domain Dq, Ai := {6 : 

Im9\, je{r) is the spherical Bessel function. 
Fix an arbitrary ^ G M.^. One can find (nonuniquely and explicitly) 9',9 E Ai, such that 
9' — 9 = ^, 9 ^ CO. For example, if = te^, t = \^\ > 0, (which can be assumed without 
loss of generality), then 9' = |e3 + ZiCi + Z2e2, 9 = — fes + ZiCi + 2:262, and the condition 

- + zl + zl = l, Zi,Z2eC, (*) 

implies 9,9' G Ai. One may find many 2:1 2:2 G C, such that Q holds and l^il 00. 
For example, take Zi = re**^, Z2 = re"*''', then r"^ sm{2ip) — sin(299) = 0, r'^[cos{2ip + 
cos(2(/?)] = 1 — ^, so cos(2(/5) = ^ — \- One can take r > \^ — |-|^/^ and find such 

that cos(2(/?) ^ ^ 



2 

In what follows we always assume 

0'^e = ^, 9',9eM, \e\^oo. (21) 
Because of (pUj) . the series 

00 

A{9',a) = J2M(^)Yi{9'), 9eM, (22) 

£=0 

converges absolutely and uniformly on compact subsets of 5*^ x A^. 

Fix positive numbers < hi < ^2 such that Dq C := {x : |x| < feo}- Note that 
the scattering solution u{x,a) for the scattering potential q{x) + C{x) can be written 
explicitly in the region |x| > 60: 

00 

m(x, a) = e*"-^' + Ae{a)Yi{a')hi{r), r := |x| > 60, (23) 
£=0 

where k = 1 and /if(r) is the spherical Hankel function which is normalized by the 
asymptotics 

hiijr) ~ — r 00. (24) 
r 

Let i^(a) G L^(S'^). Consider the problem: 

jr(z/) = min, (25) 

where 



e ^^'"^ / u{x,a)i'{a) da - I 
52 



dx. (26) 



'fei<|a;|<62 

The function u{x,a) in is defined in ()23j) . and the minimization in ()25|1 is with 
respect to all v G L^(S'^). One can prove ([TJ p. 265 ]) that 

COT} 

MT{v):=d{9)<^, 9 EM, 1^1 »1, (27) 
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Let iy{a, 9) be an arbitrary approximate solution to (j25j) in the following sense: 



(28) 



For this u define 




(29) 



where A{6',a) is defined in (j^ . 

Let ^(0 = C{x)e~'^ ''dx, where C{x) e L'^{Dq) vanishes in D'^. Let ^ hold. 

The following theorem is proved by the author in ^ p. 266]. 
Theorem. Under the above assumptions one has 



Conclusion: An algorithm is given for embedding many small particles in a domain 
Dq in such a way that the plane wave, scattered by such domain, would have a desired 
radiation pattern A{a' ^ a). 

The algorithm consists of solving the inverse scattering problem, namely, finding the 
potential q{x) + C{x), vanishing outside Dq, from the fixed-energy (fc^ = const > 0) 
scattering amplitude A{a',a). If this potential is found, then C{x) is found, and the 
small acoustically soft identical particles should be distributed in Dq with the density 
N{x) = C^^C{x), where C > / is the electrical capacitance of a perfect conductor which 
has the shape of the particle. 
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